Method for classifying multi-qubit states via the rank of coefficient matrix and 

its application to four-qubit states 
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We construct coefficient matrices of size 2 by 2 n ~ e associated with pure n-qubit states and prove 
the invariance of the ranks of the coefficient matrices under stochastic local operations and classical 
communication (SLOCC). The ranks give rise to a simple way of partitioning pure n-qubit states 
into inequivalent families and distinguishing degenerate families from one another under SLOCC. 
Moreover, the classification scheme via the ranks of coefficient matrices can be combined with other 
schemes to build a more refined classification scheme. To exemplify we classify the nine families 
of four qubits introduced by Verstraete et al. [Phys. Rev. A 65, 052112 (2002)] further into 
inequivalent subfamilies via the ranks of coefficient matrices, and as a result, we find 28 genuinely 
entangled families and all the degenerate classes can be distinguished up to permutations of the four 
qubits. We also discuss the completeness of the classification of four qubits into nine families. 



I. INTRODUCTION 

Quantum entanglement plays a crucial role in 
quantum information theory, with applications to 
quantum teleportation, quantum cryptography, and 
quantum computation The equivalence under 
stochastic local operations and classical communica- 
tion (SLOCC) induces a natural partition of quan- 
tum states. The central task of SLOCC classification 
is to classify quantum states according to a criterion 
that is invariant under SLOCC. 

SLOCC entanglement classification has been the 
subject of intensive study during the last decade 0- 
l2fj |. For three qubits, there are six SLOCC equiv- 
alence classes of which two are genuinely entangle- 
ment classes: GHZ and W 0] and four degenerate 
classes can be distinguished by the local ranks (i.e., 
ranks of single-qubit reduced density matrices ob- 
tained by tracing out all but one qubit [2] ) . For four 
or more qubits, there are infinite SLOCC classes and 
it is highly desirable to partition the infinite classes 
into a finite number of families. The key lies in find- 
ing criteria to determine which family an arbitrary 
quantum state belongs to. In a pioneering work, 
Verstraete et al. Q obtained nine SLOCC inequiv- 
alent families of four qubits using Lie group theory: 

Gabcd? Labc2i ^a^b^i L a b 3 , L a± , ^ a 2O 3 0j ) -^0503: ^^7^1 ) 

and Ln,» T n™. It is clear that, some families ob- 
tained by Verstraete et al. |3[ contain an infinite 
number of SLOCC classes and some contain both 
degenerate classes and genuinely entangled classes. 
It is of great importance to find a more refined par- 
tition of four-qubit states such that the degenerate 



classes are distinguished from the genuinely entan- 
gled families. Many other efforts have been devoted 
to the SLOCC entanglement classification of four 
qubits [5- 13|. More recently, a few attempts have 
been made toward the generalization to higher num- 



ber of qubits, including odd n qubits |17| . even n 



qubits [18|. symmetric n qubits [14H16I ] , and general 



n qubits 19|, 2j 

This paper is organized as follows. We first con- 
struct coefficient matrices of size 2 e by 2 n ~ l associ- 
ated to pure n-qubit states and prove the invariance 
of the ranks of coefficient matrices under SLOCC in 
Section II. In Section III, we present a recursive for- 
mula which allows us to easily calculate the ranks 
of coefficient matrices of n-qubit biseparable states. 
We next show that the degenerate families of gen- 
eral n qubits are inequivalent to one another under 
SLOCC in Section IV. Section V is devoted to the 
classification of four qubits via the ranks of coeffi- 
cient matrices. Section VI provides the discussion 
of the completeness of the nine families obtained by 
Verstraete et al. [3|. We finally conclude this paper 
in Section VII. 



II. THE INVARIANCE OF THE RANKS OF 
COEFFICIENT MATRICES 

Let \ip)i... n = X)i=o 1 a iV) De an "--qubit pure 
state. We associate with the state \ip)i... n a 2 l by the 
2«-f coemc i en t matrix C 1 ...^ +1 )... n (|V , )i— n) whose 
entries are the coefficients ao,a±,--- ,a2^~i of the 
state \ip)\... n arranged in ascending lexicographical 
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order. To illustrate, we list Ci...^«+i)... n (|V')i---n) 
below as: 



i n-e e Ti-t 

a o---io---o ••• a 0---ll---l 



V 



a l---10---0 



*1 - • - 1 1 • - - 1 



(1) 



In the binary form of the coefficient matrix in 
Eq. ([!]), bits I to £ and I + 1 to n are referred 
to as the row bits and column bits, respectively. If 
£ = 0, C0,i...„(|-0)i - n) reduces to the row vector 
(a ,--- ,a 2 »-i), and if £ = n, C 1 ... n ^(\ip) 1 ... n ) re- 
duces to the column vector (ao, • • • , a2"-i) T . 

Let {<7i,<72,-'' iQn} be a permutation of 
{1,2,--- ,n}. Let C qi ... qt>qe+1 ... qn (\ip) 1 ... n ) be 
the 2 e x 2 n ~ f coefficient matrix of the state \tp)i... n , 
which is constructed from the coefficient matrix 
Ci2—i,e+i— n m Eq. (H} by taking the corresponding 
permutation. Here 5i, • • • , qe are the row bits and 
qi+i,--- ,q n are the column bits. Indeed, we only 
need to specify the row bits, as the column bits 
would simply be the rest of the bits. In the sequel, 
we will omit the subscripts • • • , q n and simply 
write C qi ... qe , whenever the column bits are clear 
from the context. 

It is known that two n-qubit pure states \tp)i--- n 
and \tp')i— n are equivalent to each other under 
SLOCC if and only if there are local invertible op- 
erators Ai, A2, ■ ■ ■ , and An such that [2j 



W)x... n =Ax®A 2 ' 



•A n \ 



(2) 



In terms of coefficient matrices, it can be verified 
that the following result holds: For any two SLOCC 
equivalent n-qubit pure states \ip)i--- n and \i)')i... n , 
their coefficient matrices C qi ... qe satisfy the equation: 



C (?1 ... 3 ,(|^)l-n) = 

(A qi <g> • • • <g> A qi )C qi ... qi (\ij)) 1 ... n )(A qi+1 



invariant under SLOCC. This leads to the following 
theorem. 

Theorem 1. If two n-qubit pure states are SLOCC 
equivalent then their coefficient matrices C qi ... qi 
given above have the same rank. 

Restated in the contrapositive the theorem reads: 
If two coefficient matrices C, 



(3) 



where A\ , A2, ■ ■ ■ , and A n are the local operators in 
Eq. ([2]). Conversely, if there are local invertible op- 
erators A% , A2 , ■ ■ ■ , and A n such that Eq. ([3]) holds 



true for some C a 



{ , then \ip) i--- n and \ip')i 



equivalent under SLOCC. 

It immediately follows from Eq. ^ that the rank 
of any coefficient matrix of an n-qubit pure state is 



81— 9 



associated with 
two n-qubit pure states differ in their ranks, then 
the two states belong necessarily to different SLOCC 
classes. 

Coefficient matrices constructed above turn out 
to be closely related to reduced density matrices. 
We let pi2-n(\ip)i-n) = \ip)i-ni-n(ip\ be the den- 
sity matrix of an n-qubit pure state \tp)i... n , and we 



let p qi -.q e be the ^-qubit reduced density matrix ob- 
tained from px2— n by tracing out n — £ qubits. As 
has been previously noted for bipartite systems of 
dimensions d x d, a reduced density matrix has a 
full rank factorization in terms of the corresponding 
coefficient matrix and its conjugate transpose 1221 . 
This factorization also holds for n-qubit states [231 ] : 

Pq 1 -<u(\i>)i-n) = C qi ... qe (\ijj)i... n )Cl 1 ... qe (\ijj)i... n ), 

(4) 

where is the conjugate transpose of C . An impor- 
tant relationship between reduced density matrices 
and SLOCC polynomial invariants can be obtained 
by taking the determinants of both sides of Eq. (HJ) 
for even n and for £ — n/2, yielding: 

detpg 1 ...g n/3 (|V)i...„) = |detC gi ...< W2 (|V>)i...„)| 



Here detC, 



(5) 

(\ip)i—n) is a SLOCC polynomial 
invariant of degree 2™/ 2 for even n qubits and its ab- 
solute value can be used as an entanglement measure 
24]. Thus we have the following: 

Theorem 2. For even n-qubit pure states, the 
determinants of n/2-qubit reduced density matrices 
are the squares of the SLOCC polynomial invariants 
of degree 2™/ 2 , with the absolute values of the lat- 
ter quantifying n/2-qubit entanglement of the even 
n-qubit states after tracing out the other n/2 qubits. 

As an example, when n — 4 we have det p\2 = 
|L| 2 , detpi3 = |Af| 2 , and detp i4 = \N\ 2 , where 
L, M, and N are polynomial invariants of degree 
4 [25|. When n = 6, there are 10 three-qubit re- 
duced density matrices and 10 polynomial invariants 
of degree 8: D\,--- ,D\° For reduced density 

matrix P123 and polynomial invariant D\, we have 
det P123 = l^ejl 2 - Similar equations hold for other 
reduced density matrices and polynomial invariants 
with appropriate permutations of qubits. 



3 



Remark 1. (i). The determinants of reduced den- 
sity matrices are invariant under SLOCC. (ii). It is 
worth noting that Eq. (J5J) holds for bipartite systems 
of dimensions d x d as well [221 ] . 

As a particular case of Eq. when qi = i we 
have / Qi...„(|V>)i-ra) = Ci... n (|^)i...„)Cj... n (|^}i...„). 
By virtue of Eq. (j4} , the rank of the £-qubit reduced 
density matrix and the rank of the corresponding 
coefficient matrix are the same. In light of Theorem 
1, we have the following result. 

Corollary. The ranks of ^-qubit reduced density 
matrices obtained by tracing out n — £ qubits are 
invariant under SLOCC. 

This is particularly true for the local ranks 0]. 
Note also that any complex matrix has a singular 
value decomposition, with the number of nonzero 
singular values equal to the rank of the matrix. This 
means that the number of nonzero singular values 
of any coefficient matrix of an ?i-qubit pure state is 
invariant under SLOCC. 



III. A RECURSIVE FORMULA FOR THE 
RANKS OF iV-QUBIT BISEPARABLE 
STATES 

In principle, we can calculate the ranks of coef- 
ficient matrices for n-qubit biseparable pure states 
by direct calculations. However, in practice, this is 
rather cumbersome from the computational point of 
view, and as n becomes large, this might pose a se- 
rious problem. In order to avoid this difficulty, we 
propose a simple recursive formula for the ranks of 
n-qubit biseparable states. 

Suppose that a biseparable n-qubit pure state 
l^)i-n is of the form = \4')j 1 -j k ®\<p)j k+1 -j n 

with \4>)j 1 ---j k being a fc-qubit state and \<p)j k+1 —j n 
being an (n— fc)-qubit state. We let C qi ... qe (\ip)i... n ) 
be the coefficient matrix associated with the state 
We let C ql ... qt (\0) n ... Jk ) be the 2 s by 
2 k ~ s coefficient matrix associated with the fc-qubit 
state \4>)h-j h - Here {q*,--- ,q*} = {qi,-- - ,qe} n 
{ji , - • • ,jk} are the row bits, and by convention, 
the rest k — s bits are the column bits. Moreover, 
we let C q > i ... q > t (\cp)j k+1 ... jn ) be the 2* by 2 n - k ~ t co- 
efficient matrix associated with the (n — fc)-qubit 
state \<p)j h+1 ...j n . Here {q[, ■■■ ,q' t } = {q lr ■■■ ,qe}n 
{jk+i, ■ ■ ■ ,jn} are the row bits, and by convention, 
the rest n — k — t bits are the column bits. It can be 



verified that 

Cqi-qA\<t>)h-j k ® |p)jW"iJ 

= C qt~q%{\<t>)h-j k ) ® C q{-qi(\ l P)jl. + i~3n)- ( 6 ) 

In view of the fact that the rank of the Kronecker 
product of two matrices is the product of their ranks, 
we arrive at the following recursive formula for the 
ranks of coefficient matrices of an n-qubit bisepara- 
ble state: 

i&nk(C qi ... qe (\4>) jl ... jk ® \<p}j k+1 - jn )) 
= rank(C ?1 ....qj(|^} il ... i J)raak(Cg i ...g { (|^} iji+1 ... in )). 

(7) 

The formula above allows us to calculate recur- 
sively the ranks of coefficient matrices of n-qubit 
biseparable states in terms of the ranks of coefficient 
matrices of fc-qubit states and (n — fc)-qubit states. 
To illustrate the use of the recursive formula, we 
start with the initial values rank(CU(|(/>)yi)) = 1 and 
ra,nk(C®(\(f)) a)) = b It is known that a two-qubit 
pure state can be either of the form A-B (separable) 
or the form AB (EPR) . Using the recursive formula, 
we find rank (CU (|0) a Mb)) = ra,nk{C A (\4>) a)) x 
rank(C0(|t/9)s)) = 1. On the other hand, a direct 
calculation shows that ra,v*k(C a{W) ab)) = 2. Using 
the results obtained above, we can find the ranks of 
coefficient matrices of three-qubit pure states. Con- 
sider, for example, rank(Cc(|0)B \<p)ac)) f° r bisep- 
arable states being of the form B-AC. Using the 
recursive formula, we have rank(Cc(|</>)s \<p)ac)) — 
rank(C0(|0) B )) x rank(Cc(</?)Ac)) = 2. In a simi- 
lar fashion, we can fill in the rest of the entries in 
Table [H except those in the last row which can be 
obtained by direct calculations. Proceeding in this 
way, we can construct Tables HT1 and IIIII for the ranks 
of coefficient matrices for four and five qubits. 

Note that in Tables U and [TT| the ranks of only 
2 n_1 — 1 coefficient matrices are shown. This is due 
to the fact that interchanging two row (resp. col- 
umn) bits or exchanging the row and column bits of 
a coefficient matrix does not alter the rank of the 
matrix, since the former is equivalent to interchang- 
ing two rows (resp. columns) of the matrix and the 
latter is equivalent to transposing the matrix. Ig- 
noring C0 and C\... n which always have rank 1, this 
amounts to totally 2™ _1 — 1 potentially different co- 
efficient matrices. For example, the ranks of Cba 
and Cbc are not shown in Table Hi] since Cab and 
Cba differ by the interchange of two rows, and Cbc 
is the transpose of Cad ■ As illustrated in Tables U 
HIl and IIIII the ranks of coefficient matrices permit 
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the partitioning of the space of the pure states into 
inequivalent families under SLOCC (i.e., two states 
belong to the same family if and only if the ranks 
of coefficient matrices are all equal). In particular, 
degenerate families of three, four, and five qubits are 
inequivalent from one another under SLOCC. 



TABLE I. Ranks of coefficient matrices of three-qubit 
pure states. 



Ranks of 



Families 



Ca 



Cb 



Cc 



A-B-C 

A-BC 

B-AC 

C-AB 

ABC 



TABLE II. Ranks of coefficient matrices of four-qubit 
pure states. 



Families 


Ranks of ~ 


C B 


C c 




Cab 


Cac 


Cad 


A-B-C-D 


1 


1 


1 


1 


1 


1 


1 


A-B-CD 


1 


1 


2 


2 


1 


2 


2 


A-C-BD 


1 


2 


1 


2 


2 


1 


2 


A-D-BC 


1 


2 


2 


1 


2 


2 


1 


B-C-AD 


2 


1 


1 


2 


2 


2 


1 


B-D-AC 


2 


1 


2 


1 


2 


1 


2 


C-D-AB 


2 


2 


1 


1 


1 


2 


2 


A-BCD 


1 


2 


2 


2 


2 


2 


2 


B-ACD 


2 


1 


2 


2 


2 


2 


2 


C-ABD 


2 


2 


1 


2 


2 


2 


2 


D-ABC 


2 


2 


2 


1 


2 


2 


2 


AB-CD 


2 


2 


2 


2 


1 


4 


4 


AC-BD 


2 


2 


2 


2 


4 


1 


4 


AD-BC 


2 


2 


2 


2 


4 


4 


1 


ABCD\ 


2 


2 


2 


2 


> 2 


> 2 


> 2 



a ABCD can be further partitioned under SLOCC in 
terms of the ranks of Cab, Cac an d Cad- 



TABLE III. Ranks of coefficient matrices of five-qubit 
pure states. 



Ranks of 



Families 



C : 



i-j-k-i-m 
i-j-k-im 



i-jk-lm 



1, if a = i 

2, otherwise 



i-j-klm 

i-jklm 

ij-k£m 

ijk£m 



l u l c 

1, if a = i,j,k 1, if 0,7 = i, j,k 

2, otherwise or /3,j = £,m 

2, otherwise 

1, if P,l = 3,k 
or f3,y = £,m 

2, if /3 = i or 7 = i 
4, otherwise 

1, if a = i or j 1, if £,7 = i, j 

2, otherwise 2, otherwise 
2, if /3 = i or 7 = i 
2, 3, or 4, otherwise 

1, if /3,7 = *> j 

2, if ,8,7= M,"i 
4, otherwise 
2, 3, or 4 C 



1, if a = i 

2, otherwise 
ob 



k, £, m} is any permutation of {A, B,C, D, E}. 
' a = i, j, k, £, m. 
fit 7 = *j J; m. 



IV. DEGENERATE FAMILIES OF GENERAL 
JV QUBITS ARE SLOCC INEQUIVALENT TO 
ONE ANOTHER 



The recursive formula above further gives rise to a 
criterion for biseparability of an n-qubit pure state. 
Indeed, we note that Eq. (JTJ) holds particularly 
true for {qi,--- ,qt} = {ji,--- ,3k}- In this case, 
the coefficient matrices C 9 »... ? * and C q ^... q ' t reduce 
to a column vector and a row vector respectively, 
and therefore both of them have rank 1 . It follows 
that rank(C qi ... qe .(\(f>) qi ... qe ® \<p) qi!+1 --q n )) = L Con- 
versely, if rank(C qi ...g € (\tp)i... n )) — 1 for an n-qubit 
pure state \ip)i... n , then \ip)i... n is biseparable, being 
of the form \tp)i... n = \4>) qi ... qi ® \<p)qi +t -q n - This 
can be seen as follows. For simplicity, we assume 
qt = i with i = 1, • • • ,n. If rank(Ci 2 ...^(|^)i...„)) = 
1, then all columns of Ci%..i are proportional to 
each other and each column can be written into 
the form (aobj,aibj, ■ ■ ■ , a^-i^j) 7 Hence, \ij))\... n 
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can be written as \ip)x-- 
with \4>)i...i = Ya^q 1 a 



n = \4>)l-i <8> \<P)(t+i) 

\i)i...t and \(p)^ +1 )... n 



Sj=o ^ b)(f+i) -n- This leads to the following 
biseparability criterion for n-qubit pure states. 

Biseparability criterion for n-qubit pure states. 
For any coefficient matrix C qi ... q( associated with an 
n-qubit pure state \ip)i... ni rank(C gi ... gi (IV-')i---n)) = 
1 if and only if \ip) is biseparable, being of the form 
\ip)i... n = 10) <?!•••<« ® b) gi+ i-« n (see also [Il],[23j]). 

Invoking the fact that an n-qubit pure state is en- 
tangled if it is not full separable, we have the follow- 
ing criterion to identify n-qubit entangled (respec- 
tively, genuinely entangled) pure states: An n-qubit 
pure state is entangled (respectively, genuinely en- 
tangled) if and only if the rank of at least one of its 
coefficient matrices is (respectively, the ranks of its 
all coefficient matrices are) greater than 1. 

Note that all the above criteria can be rephrased 
in terms of the ranks of ^-qubit reduced density ma- 
trices obtained by tracing out n — £ qubits [26| or 
the number of nonzero singular values of coefficient 
matrices. 

Theorem 1 together with the biseparability crite- 
rion above yield the following theorem. 

Theorem 3. Degenerate families of general n 
qubits arc incquivalent to one another under SLOCC 
and they can be distinguished in terms of the ranks 
of coefficient matrices (or in terms of the ranks of £- 
qubit reduced density matrices obtained by tracing 
out n — £ qubits) . 

The validity of Theorem 3 can be seen as follows. 
Given an n-qubit pure state, a partition P of the n 
particles is a collection of disjoint sets in such a way 
that the particles within any one set are entangled 
and any two particles from different sets are not en- 
tangled. Suppose Fi and F 2 are two different degen- 
erate families with partitions P\ and Pi respectively. 
Without loss of generality, we assume that there ex- 
ists a set 5* such that S € Pi and S £ ?2- Then the 
states in Fx can be written in the biseparable form 
|^)s|c/?)5, where S is the set of all particles except 
those in S. According to the biseparability crite- 
rion above, rank(Cs) = 1 for states in Fx. Since the 
states in F 2 cannot be written in the above bisepa- 
rable form, rank(Cs) > 1 for states in F 2 . In light of 
Theorem 1, the two degenerate families are inequiv- 
alent to each other under SLOCC. 

In addition, we remark that degenerate families of 
general n qubits can also be distinguished from one 
another under SLOCC in terms of the ranks of £- 
qubit reduced density matrices obtained by tracing 



out n — £ qubits or the number of nonzero singular 
values of coefficient matrices. 



V. SLOCC CLASSIFICATION OF FOUR 
QUBITS VIA THE RANKS OF 
COEFFICIENT MATRICES 

Suppose that the states \ip) and \ip') of four qubits 
are SLOCC equivalent to each other, then there are 
local invertible operators Ax, A2, A3, and A4 such 
that (2j 



IV/) = Ax ® A 2 <8> ^3 ® A 4 



(8) 



For a four-qubit state = Xaio ai N)' we consider 
three coefficient matrices Cab, Cac, and Cad as 
follows: 



Cai 



Cac 



C a 



D 





«() 


ax 


« 2 


03 \ 




a 4 


a 5 


0,5 


a 7 




as 


« 9 


aio 


ax\ 


V 


a\2 


ai 3 


ai4 


ax5 J 




«o 


ax 


a 4 


as \ 




« 2 


as 


«6 


a 7 




a.% 


09 


ax2 


ax3 


V 


aw 


ax\ 


ai4 


axs J 


/ 


a 


«4 


"2 


a& \ 




ax 


O5 


"3 


a 7 




as 


a\2 


aio 


a i4 


V 


09 


ai3 


an 


ais / 



(9) 



(10) 



(11) 



The coefficient matrices above satisfy the following 
equations: 

C AB {W)) =Ax® A 2 C A B(\i>))(A 3 ® Aa) t , (12) 
C AC {W)) = Ax®A 3 C A c(m(A2® A 4 f,(13) 
Cad{W)) =Ax® A4C AD (m)(A 3 ® .4 2 ) T (14) 



It follows from Eqs. (JT2J) - (JT4J) that if two four- 
qubit states are SLOCC equivalent then their coef- 
ficient matrices Cab (and also Cac and Cad) have 
the same rank. Conversely, if one of the coefficient 
matrices Cab, Cac, and Cad differ in the ranks, 
then the two four-qubit states are SLOCC inequiv- 
alent. Let family F^J 3 be the set of all four-qubit 
states with the same rank tab of the coefficient ma- 
trix Cab- Here tab ranges over the values 1, 2, 
3, and 4. Clearly, each one of the nine families 
introduced by Verstraete et al. [3] can be further 
divided into four SLOCC inequivalent subfamilies 
corresponding to the four possible values of tab ■ In 
a similar manner, we can define the families F^j^? 



G 



and F^ A A ° . One can obtain a more refined parti- variant theory. Let L' ab be defined by 
tion by further dividing the families FrA 3 , F?/£ 



4C 

TAB 1 * r A C ' 

and Ff/° into subfamilies F^GacCad _ f Cab n 
F?S°nF%%>. Clearly, the subfamilies F^ r c A ^A D 
and f c abCacC A d are SLQCC inequivalent when 

V AB r AC V AD 

TABTAcrAD ^ r' AB r' Ac r' AD . 

We now further partition the nine families intro- 
duced by Verstraete et al. 0] into SLOCC inequiv- 
alent subfamilies via the rank of coefficient matrix. 
For convenience, we rewrite the families G a bcd and 
L a bc 2 as: 

G abcd = a(\0) + |15)) + 0(13) + |12» + 7 (|5) + |10)) 
+5(|6) + |9)), (15) 

L abC2 = a'(\0) + |15» + /3'(|3) + |12)) + 7 '(|5) + |10>) 
+|6). (16) 

In Table HVl we show the subfamilies F^ AB , 

<——• TAB ' 

Ff A A c c , and F^/° of G abcd . As illustrated in Table 
fVl Gabcd can be further partitioned into nine gen- 
uinely entangled subfamilies and three biseparable 
subfamilies (marked with "*") via tab, tac, and 
tad (subfamilies not listed in the table are empty). 
For simplicity, the detailed descriptions of the sub- 
families are not shown as they can be easily obtained 
by taking the intersections of the corresponding de- 
scriptions in Table El Tables EH and EED illustrate 
the partitions of the other eight families introduced 
by Verstraete et al. into inequivalent subfamilies. In 
total, we find 28 genuinely entangled subfamilies and 
all the degenerate classes can be distinguished up to 
permutations of the four qubits (i.e., A-B-C-D, A- 
B-CD,AB-CD, \0) A \W) BCD, and 1 0) A | GHZ) BCD)- 



L' ab3 = a(|0000) + |1111)) 



a + b 



(|0101) + |1010)) 



a — b 



-(|0110) + |1001» 



+ -W|0001) + |0010) - |0111) - |1011))(18) 
v2 

that is, L' abs is obtained by replacing the two "+" 
signs of the last two terms in the formula of L a03 
by "-" signs 0|. It is claimed that there is a perfect 
correspondence between the nine families obtained 
by Verstraete et al. (with L ab3 replaced by L' ab ) 
and the nine families obtained by Chterental et al. 
[6j|. Note that the formula of L' h has also been 



adopted in Ref. [ll|. Since both Verstraete et al. 



and Chterental et al. claimed that the nine families 
obtained in their work are inequivalent to each other, 
a detailed study of the relation between L ab , d and 
L' ab can provide insights into the completeness of 
their classifications. 



A. L a b 3 (a = 0) is SLOCC equivalent to L' aba (a = 0) 

It is readily verified that the following equation 
holds between L' ab3 (a = 0) and L a b 3 (a = 0): 

L 'ab 3 = °) = 1 ® 1 ® ia z ® iv z L ab3 (a = 0), (19) 

where / is the identity and a z = diagjl, —1}. 

It follows from Eq. (fT9| that L a b 3 (a — 0) and 
L' ab (a = 0) are SLOCC equivalent. In particu- 
lar, setting 6 = yields that the states 4=(|0001) + 
|0010) - |0111) - 1 1011)) and -^(|0001) + |0010) + 
|0111) + |1011» are equivalent under SLOCC. 



VI. DISCUSSION OF THE COMPLETENESS 
OF THE NINE FAMILIES OBTAINED BY 
VERSTRAETE ET AL. 



B. L' ab3 (a ^ 0) [respectively, L aH (a / 0)] is 
SLOCC inequivalent to L a t, 3 (respectively, L' ab ) 



The family L a b 3 in Ref. [3j was defined as 
L ab3 = a(|0000) + |1111» + lt-(|0101) + |1010)) 



a — b 



(|0110) + |1001)) 



+-=(|0001> + |0010) + |0111) + |1011)X17) 
v2 



In later work, Chterental et al. [3j obtained nine 
SLOCC inequivalent families of four qubits using in- 



We first show that the family L' ab3 (a ^ 0) is 
SLOCC inequivalent to the family L ab3 . In Table 
IVIIII we show the partition of L' ab3 into SLOCC in- 
equivalent subfamilies via tab, tag-, and tad- Con- 
sulting Tables I VIII and IVIIII and using the fact that 
the subfamilies with different ranks of coefficient ma- 
trices are SLOCC inequivalent to each other, it suf- 
fices to consider the following six cases. 

Case 1. L' ab (a = 6^0) is SLOCC inequivalent 
to L ab3 {b = -3a ^ 0). 
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TABLE IV. The subfamilies Fr/ B B , F°/ C c , and Fr A A £ of G abcd . 



Subfamily 


Description 




f uab 


a = /3 = 0& 7 = ±<5^0| 


a = ±/3^0&i-y = S 




q = /? = 0& 7^±(5|7 = 


<5 = 0& a/±/3|a 






= ±(5 # & a # ±/3 


p^AB 


a ± ±/3 & 7 / ±5 




pC AC 


a = 7 = 0& /? = ±<5/0 


a = ±7^0&/3 = <5 


JpCAC 
r 2 


a = 7 = 0& /3^±(5|/3 = 


5 = 0& a/±7|a 


TpCAC 
r 3 


a = ±7^0&/3/±<5|/3 


= ±(5/0&a^±7 


t?Cac 
r 4 


a # ±7 & /3 # ±(5 




j?Cad 


a = (5 = 0& /3 = ±7/0 


a = ±5/0&/3 = 7 


pCAD 


a = 5 = 0& /3/±7|/? = 


7 = 0& a^±<5|a 


P^AD 


q = ±5^0&/3/±7|/3 


= ±7 / & a ± ±8 


P^AD 


a / ±5 & /3 / ±7 





= 

±/3/0&7 = ±5#0 



= 

± 7 /0&/3 = ±<5/0 



= 

±<5^0&/3 = ±7/0 



TABLE V. SLOCC classification of Gated via tvlb, tac, 
and tad- The subfamilies marked with "*" are bisepa- 
rable. 



TAB VAC TAD 


Subfamily description 


222 




n F% AC n F% AD 


244 


P&AB 


n f 4 Cac n F° AD 


333 


j?Cab 


n F° AC n F° AD 


344 


T?C AB 


n Ff- 40 n f1 ad 


424 


P&AB 


n Ff- 40 n f 4 Cad 


434 


j?Cab 


n F° AC n f 4 Cad 


442 


t?Cab 
r 4 


n f 4 Cac n f 2 Cad 


443 


P^AB 


n f 4 Cac n F° AD 


444 


pC AB 
r 4 


n f 4 Cac n F° AD 


144* 


pCAB 


(i.e., AB-CD) 


414* 


pCAC 


(i.e., AC-BD) 


441* 


P&AD 


(i.e., AD-BC) 



In this case, we can resort to -D™, a degree 6 
polynomial invariant of four qubits [25| (see the Ap- 
pendix for the expression of D xy ). Indeed, it can 
be verified that if \ip) and are any two SLOCC 
equivalent states, that is, they satisfy Eq. ([2]), then 
the following equation holds: 

D xy (W)) = D xv {\ip))[U A t=l detA l f. (20) 

It follows from Eq. (JSC)) that for any two SLOCC 
equivalent states and \ip'), either D xy (\tp')) and 
D xy (jip)) both vanish or neither vanishes. 



A direct calculation shows that 

D xy = -^(a-bf(a + b) 3 (21) 

for both L a f, 3 and L' ab . The desired result then fol- 
lows by noting that D xy = 16a 6 7^ for L a b 3 (b — 
—3a 7^ 0) whereas D xy = for L' ab:t (a = b 7^ 0). 

Case 2. L' ab3 (a = — b 7^ 0) is SLOCC inequivalent 
to L ab:i (b = 3a 7^ 0). 

This case can be dealt with similarly as case 1 by 
noting that D xy = 16a 6 7^ for L a b 3 (b — 3a 7^ 0) 
whereas .D^ = for L' ab3 (a = — b 7^ 0). 

Case 3. 1/ h (6 = —3a 7^ 0) is SLOCC inequiva- 
lent to L a b 3 {b = —3a 7^ 0). 

In this case, the semi-invariants defined in Ref. [Tj 
turn out to be useful. More specifically, for any four- 
qubit state = J2i=o Ci\i), the semi-invariants F± 
and F2 are defined in Ref. [7j as 

Fi(ip) = (c c 7 - c 2 c 5 + cic 6 - C3C4) 2 

-4(c 2 c 4 - c c 6 )(c 3 c 5 - C1C7), (22) 

^(VO = (C8C15 - C11C12 + C9C14 - C10C13) 2 

-4(ciici 3 - c 9 ci 5 )(cioCi2 - c 8 ci 4 ).(23) 

Let |0) be any four-qubit state SLOCC equivalent 
to L a b 3 [i.e., they satisfy Eq. @]. Let 

A, = ( ai a2 ). (24) 
\a 3 a 4 J 
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TABLE VI. SLOCC classification of L a bc 2 via tab, tac, and tad- The subfamilies marked with "*" are biseparable. 



tab tac tad Subfamily description 



233 


a! 


= P' = 


-- o & i ± o 




244 


a' 


= ±p' =/= o & y = o 




323 


a' 


— 1 — 


= o & p' ^ 




332 


a' 


& 


^ = y = o 




333 




= ±P' 


= ±7 V 




344 


i 


= & 


q'/3' ^ & a' # ±/3' | 7' ^ & a' 


= ±/3' / & a V ±7' 


424 


P' 


= & 


a' = ±7 V 




434 


P' 


= & 


a'7' ^ & a' # ±7' | P' & a' 


= ±7' ^ & a' ^ ±/8' 


442 


a' 


= & 


p' = ±y ^ o 




443 


a' 


= & 


/3' / ±7' & /3' 7 ' / | a V & /3' 


= ±7' / & a' / ±/3' 


444 


i 


/ & 


a V ±/9' & P V & a V ±7' & a' / & /3 V ±7' 


111* 


a! 


=p' = 


= 7' = 0, (i.e., A-B-C-D) 





TABLE VII. SLOCC classifications of L ab3 , L a2 t 2 , L ai , L a2 o 3S)I , £o 5e3 , Lq 7@1 , and Lo 3eI o 3e i via tab, tac, and 
tad- The subfamilies marked with "*" are biseparable. 



Family tab tac tad Subfamily description 



Family 



tab tac tad Subfamily description 



Lao 



La 4 



333 

424 

434 

212* 

323 

434 

333 

222* 

333 



ab = & a ± b 

a = ±6 / 

ab / & a / ±b 

a = b = (i.e., A-C-BD) 

L aA (a = 0) 

L a4 (a^0) 

L «20 3eI (a / 0) 

a = (i.e., \Q)a\W)bcd) 

Lo vml 



Lab 



222 


a = b — (i.e., 


344 


ab = & a / 6 


424 


a = b / 


434 


& = -3a / 


442 


a = -b=£Q 


443 


b = 3a / 


444 


ab / & b ± ± 


333 




222* 


|0}a|GHZ> S cd 



TABLE VIII. SLOCC classification of L' aH via tab, tac, 
and tad- 



TAB TAC TAD 



222 
344 
424 
434 
442 
443 
444 



Subfamily description 



a = b = (i.e., |W)asod) 

ab = & a / 6 



a = / I b = -3a / 


a = -o/0|fc = 3a/0 
a6/0&6/±a&6 / ±3a 



A tedious but straightforward calculation yields 
1 



r 4 



F 2 (^) = -(a 2 -6 2 )^ 



\\ det At 

1=2 
4 

det Ai 



(25) 
(26) 



In view of Eqs. (|23|) and (|2l)f and the fact that „4i 
is invertible, it follows at once that if \(f>) is SLOCC 
equivalent to L a h 3 (a ^ ±&), then the following equa- 
tion holds: 



l^i Wl + \F 2 (<t>)\ ? 0. 



(27) 



Let |iy9) be any state SLOCC equivalent to L' Qb3 
[i.e., they satisfy Eq. ||5J)]. Again, a tedious but 
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straightforward calculation yields 



Fib) 



-1 



2^2 



iV2(3a 2 + b 2 )ai + 8a(a 2 - b 2 )a 2 ) 



JJdetA 



*2(¥>) 



L i=2 
-1 



(28) 



2^2 

r 4 



ia|(-iV2(3a 2 + o 2 )a 3 + 8a(a 2 - 6 2 )a 4 ) 



j J det Ai 



When a (a 2 — b 2 ) 7^ 0, consider the operator 



\ a 4 



(29) 



(30) 



where aia 4 7^ 0. Clearly, A\ is invertible. In view 
of Eqs. ([25 ]) - ([50 ]) . it follows that there exists a state 
\ip*) equivalent to L' ab3 (a(a 2 — b 2 ) 7^ 0) under local 
invertible operators A\, A%, A3, and A4, such that 



|iq(^*)| + |F 2 (^*)|=0. 



(31) 



From Eqs. (27J) and ([3TJ), is SLOCC in- 

equivalent to the state L a b 3 (a 7^ ±6). Therefore, 
L 'ab 3 ( a ( a2 - b2 ) ^ °) is SLOCC inequivalent to 
L a b 3 {oL 7^ ±6). In particular, L' ab (b = —3a 7^ 0) 
is SLOCC inequivalent to L a b 3 (b = — 3a 7^ 0). 



Case 4. L' , 



3a 7^ 0) is SLOCC inequivalent 



to L ab3 (b = 3a 7^ 0). 

This case can be treated analogously to case 3. 

Case 5. L' ab3 (a ^ & b = 0) is SLOCC inequiva- 
Jent to L a f, 3 (a5 = & a 7^ 6) . 

In Ref. [HI, we proved that L a & 3 (a = & 6 7^ 0) 
and L a b 3 (a 7^ & b — 0) are SLOCC inequivalent. 
A proof analogous to that of Ref. [l(| shows that 
L' ab3 (a = & b 7^ 0) and L' Q&3 (a 7^ & & = 0) are 
SLOCC inequivalent. Using the fact that L a f, 3 (a — 
& 6 7^ 0) is SLOCC equivalent to L' abi (a = & 6 7^ 
0) [see Eq. (dSJ] yields that L^ 3 (a 7^ & 6 = 0) 
is SLOCC inequivalent to L a b 3 (a = & 6 7^ 0). 
Furthermore, an argument analogous to case 3 shows 
that L' ab (a 7^ & b = 0) is inequivalent to Z (, 3 (a 7^ 
0&6 = 0). 

Indeed, we can further conclude that L a b 3 (a = 
0) and L a b 3 (a 7^ 0) are SLOCC inequivalent and 
L' ab3 (a — 0) and L' ab (a 7^ 0) are SLOCC inequiva- 
lent. 



Case 6. L' ab3 (ab ^0ka^±bkb^ ±3a) is 
SLOCC inequivalent to L a b 3 (ab 7^ & a 7^ ±6 & 7^ 
±3a). 

This case can be treated analogously to case 3. 

As a consequence, L' ab3 (a 7^ 0) is SLOCC inequiv- 
alent to L a b 3 . An analogous argument shows that 
L a b 3 (a 7^ 0) is SLOCC inequivalent to L' ab . 



C. The relation between L' ab3 and L a b 3 under 
per mutat ions 

Let 1 7) be the state of the subfamily L' ab3 (a 7^ 
& b — 0), 1 17) be the state of the subfamily 
L' ab3 (b = 3a 7^ 0), be the state of the sub- 
family L' ab (b = —3a 7^ 0), and \v) be the state 
of the subfamily L' abs (ab ^0ka^±bkb^ 
±3a). We argue that the above four subfamilies 
are SLOCC inequivalent to L a b 3 under any permu- 
tation of qubits. This can be seen as follows. Let 
be the transposition of qubits i and j. A te- 
dious calculation shows that the permutations giv- 
ing rise to different (7) are k,\ = 7, K2 = (1,3), 
k 3 = (1,4), k 4 = (1,2)(1,3), k 5 = (1,2)(1,4), and 
K6 = (1,4)(1,2)(1,3). Similarly, the permutations 
giving rise to different \rj), I'd), and \v) are 7Ti = 7, 

7T 2 = (1,2), 7T 3 = (1,3), 7T 4 = (1,4), TT 5 = (1,3)(1,2), 

7T 6 = (1,4)(1,2), vr 7 = (1,2)(1,3), vr 8 = (1,2)(1,4), 
7T 9 = (1,2)(1,3)(1,2), vr 10 = (1,2)(1,4)(1,2), tt u = 
(1,4)(1,2)(1,3), and tt 12 = (1, 4)(1, 2)(1, 3)(1, 2). 
The result that Ki\j)(i — 1, • • • ,6), tt^ 1 77} , 7Tj|$), and 
Kjl^U = 1, • ■ ■ 1 12) are all SLOCC inequivalent to 
L a b 3 then follows by calculating the ranks Tab-, tac, 
and tad of ^7), Tj|?7)j tt j | i?) and tTj | , and using 
an argument analogous to that of case 3 in the pre- 
vious section. 

Remark 2. By using Tables IVIII and IVIII1 one can 
verify that (l,4)L' o6s (o = b ^ 0) is SLOCC equiva- 
lent to L a b 3 (a = & b 7^ 0) under the invertible local 
operator a x ®a z ®il®a y , and (1, 3)L' ab (a = — b 7^ 0) 
is SLOCC equivalent to L a b 3 (a = & b 7^ 0) under 
the invertible local operator a x ® cr z ® cr y ® il. 



D. I/^f, 3 (a 7^ 0) is SLOCC inequivalent to the 
other eight families by Verstraete et al. 

Here we show that L' ab3 (a 7^ 0) is not only SLOCC 
inequivalent to L a & 3 but also SLOCC inequivalent 
to the other eight families by Verstraete et al. For 
simplicity, we only show that L' ab (a = — b =/= 0) 
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is SLOCC inequivalent to the other eight families 
obtained by Verstraete et al. From Table IVIII1 
tab? Ac? ad = 443 f or L' ba (a = -b ^ 0). Consult- 
ing Tables[V]|VTj andlVII| and using the fact that the 
subfamilies with different ranks of coefficient matri- 
ces are SLOCC inequivalent to each other, it suffices 
to show that L' abs (a = — b ^ 0) is SLOCC inequiv- 
alent to the subfamilies with r abt 'ac t 'ad — 443 of 
Gabcd and L a bc 2 ■ 

To show that L' ab3 {a =-6^0) is SLOCC in- 
equivalent to the subfamily with r abt Ac r ad = 443 
of Gabcd, we use the degree 6 polynomial invariant 
D xy given in Eq. (f20|l . It is readily seen from Eq. 



jH} that D xy = for L' at>3 (a = -b ^ 0). A simple 
calculation shows that 

D xy = { a p- 1 5)(al3 + 1 8){a 2 +f3 2 - 1 2 - 5 2 ) (32) 

for Gabcd [as defined in Eq. ([To]) ]. It is readily seen 
from Eq. (I3"2"j) that D xy 7^ for the subfamily with 
t abt ac t ad — 443 of Gabcd and then the desired 
result follows. 

Next we show that L' aba (a = -6^0) is SLOCC 
inequivalent to the subfamily with tab^actad = 
443 of L a bc 2 [ as defined in Eq. ((To]) ]. A calculation 
shows that 



D xy = {a'li') 2 ( a ' 2 -i 2 +f5' 2 ) 



(33) 



for L a b C2 . From Table [VH we distinguish the follow- 
ing two cases. 

Case 1. a' ^ k ? = ±7' ^ k a' ^ ±f3' . 
In this case D xy 7^ and then the desired result 
follows. 

Case 2. a' = k ft ^ ±7' k /3' 7 ' ^ 0. 
In this case D xy ~ 0. We can resort to the semi- 



invariants given in Eqs. (|2"2"|) and (1231) . Let \ip) be 
any state SLOCC equivalent to L' ab3 (a = — b 7^ 0) 
with A given by Eq. (|24|) . A tedious but straight- 
forward calculation yields 



F 1 (|^)) = -2a 2 af 



F 2 (\<p)) = -2c?a% 



Y\ dot A 

t=2 
4 

n dct ^ 



n 2 



(34) 
(35) 



In view of Eqs. ([3~4l and ([33)) and the fact that A 
is invertible, it follows at once that if \(f) is SLOCC 
equivalent to L' ab3 {a = — b 7^ 0), then the following 
equation holds: 



(36) 



The desired result then follows by noting that i*i = 
F 2 = for L abC2 with a' = k /?' 7^ ±7' & /3' 7 ' ^ 0. 

As a consequence, L' ab3 (a = — b 7^ 0) is SLOCC 
inequivalent to the nine families obtained by Ver- 
straete et al. Q. 

The discussion suggests that the partition in Ref. 
Q is incomplete. For completeness, one may add 
the family L' abs to the family L a b 3 in Ref. [3j. An 
analogous argument shows that the partition in Ref. 
6j is incomplete as well, and for completeness, one 
may add the family L a b 3 to the family 6 in Ref. Q . 



VII. CONCLUSION 

We have recast the necessary and sufficient con- 
dition for two n-qubit states to be equivalent under 
SLOCC into an equivalent form in terms of the co- 
efficient matrices associated with the states. As a 
direct consequence of the new necessary and suffi- 
cient condition, we have showed that the rank of the 
coefficient matrix as well as the rank of the £-qubit 
reduced density matrix is invariant under SLOCC. 
We have also presented a recursive formula for the 
calculation of the rank of coefficient matrix of an n- 
qubit biseparable state. The recursive formula fur- 
ther gives rise to a biseparability criterion in terms 
of the rank of coefficient matrix to determine if an 
arbitrary n-qubit pure state is biseparable. The in- 
variance of the rank of coefficient matrix together 
with the biseparability criterion reveals that all the 
degenerate families of general n qubits are inequiv- 
alent under SLOCC. 

We have then classified four-qubit states under 
SLOCC via the ranks of coefficient matrices and the 
nine families introduced by Verstraete et al. were 
further partitioned into inequivalent subfamilies. In 
particular, we have found 28 genuinely entangled 
families and all the degenerate classes can be dis- 
tinguished up to permutations of the four qubits. 
We have performed a detailed study of the rela- 
tion between the family L a b 3 and the family L' at)3 
with corrections to the signs of the last two terms 
in the formula of L a b 3 via the ranks of coefficient 
matrices. By using a degree 6 polynomial invari- 
ant and two semi-invariants of four qubits, we have 
found that L' ab3 (a = 0) is SLOCC equivalent to 
L' ab (a = 0) whereas L' ab (a 7^ 0) is SLOCC inequiv- 
alent to L a03 {a 7^ 0). We have also demonstrated 
that L' ab3 (a ^ k 6 = 0), L' ab3 {b = ±3a 7^ 0), and 
L' ab3 {ab ^0ka^±bkb^ ±3o) are SLOCC in- 
equivalent to L a b 3 under any permutation of qubits, 
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whereas L' ab3 (a = ±b ^ 0) are SLOCC equivalent 
to L a b 3 (a — & b 7^ 0) under some permutations. 
This suggests that the partition of four-qubit states 
into the nine families by Verstraete et al. is incom- 
plete, and for completeness, one may simply add the 
family L' ab to the family L all3 . 
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APPENDIX 



Following [25j, D xy can be constructed as 





dn 


d\2 


du 


D X y — 




d 2 2 


d 23 




d 3 i 


d 3 2 


d33 



(37) 



where the entries of D xy are given by 



di2 
du 
d 2 i 
d 22 

d 2 3 

d 3 i 
d 32 

d33 



V 

a a 3 - aia 2 , 

aoa7 — ciiae — a 2 a 5 + 0304, 
a 4 o 7 - a 5 a 6 , 

a an - aia w - a 2 a 9 + a 3 a 8 , 
ao a i5 — ai<ii4 — a 2 a\3 + a 3 ai2 
+a 4 an - a 5 a w ~ a 6 a 9 + a 7 a 8 , 
a 4 a 15 - a 5 a 14 - a 6 a 13 + a 7 a 12 , 

08011 — 09^10) 

08«15 — a 9 ai4 — a 10 ai3 + anai2, 
012015 — 013014. 



(38) 
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